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Verlinde's proposal on the entropic origin of gravity is based strongly on the assumption that 
the equipartition law of energy holds on the holographic screen induced by the mass distribution 
of the system. However, from the theory of statistical mechanics we know that the equipartition 
law of energy does not hold in the limit of very low temperature. Inspired by the Debye model 
for the equipartition law of energy in statistical thermodynamics and adopting the viewpoint that 
gravitational systems can be regarded as a thermodynamical system, we modify Einstein field equa- 
' tions. We also perform the study for Poisson equation and modified Newtonian dynamics (MOND). 

Interestingly enough, we find that the origin of the MOND theory can be understood from Debye 
entropic gravity perspective. Thus our study may fill in the gap existing in the literature under- 
standing the theoretical origin of MOND theory. In the limit of high temperature our results reduce 
to their respective standard gravitational equations. 
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Thermodynamics of black holes reveals that geometrical quantities such as horizon area and surface gravity are 



Y\ ' related to the thermodynamic quantities such as entropy and temperature. The first law of black hole thermodynamics 

. ,— i \ implies that the entropy and the temperature together with the energy (mass) of the black hole satisfy dE = TdS [l| . In 

!*^| ■ 1995 Jacobson [2[ put forwarded a new step and suggested that the hyperbolic second order partial differential Einstein 

c* *', equation for the spacetime metric has a predisposition to thermodynamic behavior. He disclosed that the Einstein 

O . field equation is just an equation of state for the spacetime and in particular it can be derived from the proportionality 

i— i' of entropy and the horizon area together with the fundamental relation 8Q — TdS. Following Jacobson, however, 



several recent investigations have shown that there is indeed a deeper connection between gravitational dynamics and 
horizon thermodynamics. The deep connection between horizon thermodynamics and gravitational dynamics, help to 
understand why the field equations should encode information about horizon thermodynamics. These results prompt 
£f*l • people to take a statistical physics point of view on gravity. 

A next great step put forwarded by Verlinde [3j who claimed that the laws of gravity are not fundamental and in 
. particular they emerge as an entropic force caused by the changes in the information associated with the positions of 
' material bodies. According to Verlinde proposal when a test particle with mass m approaches a holographic screen 
: from a distance Ax, the change of entropy on the holographic screen is 

^ : 

^ . AS = 2tt— Az, (1) 

where we have set fee = c = 1 for simplicity, through this paper. The entropic force can arise in the direction of 
increasing entropy and is proportional to the temperature, 

F = T^. (2) 

Ax y ' 

Verlinde's derivation of Newton's law of gravitation at the very least offers a strong analogy with a well understood 
statistical mechanism. Therefore, this derivation opens a new window to understand gravity from the first principles. 
The study on the entropic force has raised a lot of attention recently (see fil-fll| and references therein). 

Verlinde's proposal on the entropic origin of gravity is based strongly on the assumption that the equipartition 
law of energy holds on the holographic screen induced by the mass distribution of the system, namely, E — \NT. 
However, from the theory of statistical mechanics we know that the equipartition law of energy does not hold in the 
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limit of very low temperature. By low temperature, we mean that the temperature of the system is much smaller than 
Debye temperature, i.e. T <C T D . It was demonstrated that the Debye model is very successful in interpreting the 
physics at the very low temperature. Hence, it is expected that the equipartition law of energy for the gravitational 
systems should be modified in the limit of very low temperature (or very weak gravitational field) . 

It is important to note that Verlinde got the Newton's law of gravitation, Einstein equations and Poisson equation 
with the assumption that each bit on holographic screen is free of interaction. It should be more general that the bits 
on holographic screen interact each others. In such case, one could anticipate that the Newton's law of gravitation, 
Einstein equation and Poisson equation must be modified. For example, Gao [l2| studied three dimensional Debye 
model and modified the entropic force and henece Friedmann equations. Such modification can interpret the current 
acceleration of the universe without invoking any kind of dark energy [l2| ■ In this paper we use the Debye model to 
modify the entropic gravity. We find that this modified entropic force affects on the law of gravitations and modify 
them accordingly. 

This paper is structured as follows. In the next section we derive Einstein field equations from Debye entropic 
gravity. The theoretical origin of MOND theory is discussed in the framework of Debye entropic gravity in section 
III. Sec. IV is devoted to the derivation of the Poisson equation from Debye entropic force scenario. We finish our 
paper with conclusions which appear in Sec. V. 



II. EINSTEIN EQUATIONS FROM DEBYE ENTROPIC GRAVITY 

Following Verlinde's scenario, gravity may have a statistical thermodynamics origin. Thus, any modification of 
statistical mechanics should modify the laws of gravity accordingly. In this section we use the modified equipartition 
law of energy to obtain the modified Einstein equations. 

We consider a system that its boundary is not infinitely extended and forms a closed surface. We can take the 
boundary as a storage device for information, i.e. a holographic screen. Assuming that the holographic principle 
holds, the maximal storage space, or total number of bits N, is proportional to the area A, 

Suppose there is a total energy E present in the system. Let us now just make the simple assumption that the 
energy is divided evenly over the bits. Each bit on the holographic screen has one dimensional degree of freedom, 
hence we can use the one dimensional equipartition law of energy. The equipartition law of energy which is valid in 
all range of temperatures is 

E = \nTD{x), (4) 
where T is the temperature of the screen and D(x) is the one dimensional Debye function defined as 

D {x) ^-f Q ^- ld y, (5) 

and x is related to the temperature 

* = (6) 

where To is the Debye temperature. Using the equivalence between mass and energy, E — M, as well as Eq. ([3]), we 
can rewrite Eq. ((4]) in a more general form, 

where the integration is over the holographic screen. For temperature, we use the Unruh temperature formula on the 
holographic screen, 

ha , . 

T= ^' (8) 

where a denotes the acceleration. The acceleration has relation with the Newton's potential and in general relativity 
it may be written as 

-b 



a b = - VV, (9) 
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where <j> is the natural generalization of Newton's potential in general relativity and for it we have [l3| , 

0=\Ln(-CHa), (10) 

where £° is a global time like Killing vector. The exponent represents the redshift factor that relates the local time 
coordinate to that at a reference point with <fi = 0, which we will take to be at infinity. We choose the holographic 
screen S as a closed equipotential surface or in other words, a closed surface of constant redshit <j>. Therefore Eq. (jHJ 
may be written as Q 

T=A e ^v a V a 0, (11) 

where iV a is the unit outward pointing vector that is normal to the equipotential holographic screen S and time 
like Killing vector £ b . We inserted a redshift factor e^, because the temperature T is measured with respect to the 
reference point at infinity. Because N a is normal to the equipotential holographic screen, for it we have 

W (12) 



Therefore we can rewrite Eq. (jlip as 



Substituting Eq. ((TTJ» in Eq. flT}, we get 



T= ^(WV^) 1 / 2 . (13) 



M = — !— * e' t 'N a V a ^D(x)dA. (14) 
47rG Js 

Following the same logic of [l3| , we can obtain 

M = -^j>V a l*D{x)dS ab , (15) 
where d5 n b is the two-surface element [l4| • On the other hand, according to the Stokes theorem, we have [l4[ 

B ab dS ab = 2 f V b B ab dZ a , (16) 

S JY. 

where B ab is an antisymmetric tensor field and S is the two dimensional boundary of the hypersurface S. dT, a is a 
directed surface element on S and for it we have 

cE Q = en a dT,, (17) 

where n a is the unit normal of the hypersurface E and e is equal to -1 or 1 if the hypersurface is spacelike or timelike, 
respectively. Now we apply the Stokes theorem (|16l) for Eq. (TT51 and get 

M = j V h [V a i h D{x)]d^ a 

= J[D(x)V b V a e + V a £ b V b D(x)]dE a 

= J[-D(x)V b V b r + V a e^ b D(x)]dJ: a , (18) 

where in the last step we have used the Killing equation, 

V a £ fc + V fc r = 0. (19) 

Now we use the relation [l]| 

V a V Q £" = -R b J a , (20) 
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which is implied by the Killing equation for £ a , and get 

M = JjR ab ^ b D(x) + V a £ c V c D(x)}dZ a 

= J [R a ^ h D{x) +e-^{~e^a^cD{x)]dY. a 

= -^qJ^D{x) - e- 2 %V a Z c V c D(x)]n a Z b dE, (21) 

where in the second line we have used Eq. (|TU| . In the last line we have used dS a = — n a dE, because the hypersurface 
£ is spacelike. 

On the other hand, M can be expressed as an integral over the enclosed volume of certain components of stress 
energy tensor Tab UM, 

M = 2 j {Tab - ^Tg ab )n a edZ. (22) 

Equating Eqs. flU) and ([22]), we find 



D(x)R ab - e- 2 *tbVat; c VcD{x) = ^G(Tb - \r 9 ab)- (23) 

The above equation is the modified Einstein equations resulting from considering the Debye correction to the equipar- 
tition law of energy in the framework of entropic gravity scenario. This equation is now valid for all range of 
temperature, since we have assumed the general equipartition law of energy. Therefore, we see that in Verlinde's 
approach, any modification of first principles such as equipartition law of energy will modify the gravitational field 
equations. The question whether the modified term in Einstein equation can be detectable practically or not needs 
more investigations in the future. One needs to first specify the Debye function D(x) and then try to solve the field 
equations (|2"3"|) . The resulting solutions should be checked with experiments or observations. It is clear that the 
correction term only plays role in very low temperature, in which the curvature of spacetime tends to zero and it 
becomes flat. 

It is instructive to examine the modified Einstein equations in the high temperatures limit. According to the Unruh 
temperature formula we have 

9 = fr, (24, 

where g is the norm of the gravitational acceleration. Therefore, the strength of the gravitational field is proportional 
to the temperature. Also, we can define the Debye acceleration relating to the Debye temperature as 

9D = ^T D . (25) 
n 

Therefore, if the temperature is larger than the Debye temperature, i.e. T > To, then the norm of the gravitational 
acceleration is larger than the Debye acceleration, i.e. g > go- In other words, the limit of high temperatures 
compared to the Debye temperature, is corresponding to the strong gravitational fields. In this case we have T 3> Tq, 
thus for x and y in the definition of the Debye function ([3]), we have i < 1 and consequently !/ C 1. Therefore we 
can use the approximation e v ~ 1 + y in the integral of Eq. ([5]) and as a result, the one dimensional Debye function 
reduces to 

1 r 

D(x) w - / dy = l. (26) 
x Jo 

Substituting this result (D(x) — 1) in the modified Einstein equations (|2"3"T) . leads to 

Rab = SirGiTb - \Tg ab ). (27) 



Therefore, in the temperatures extremely larger than the Debye temperature (very strong gravitational fields), one 
obtains the standard Einstein field equations as expected. 
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III. MOND THEORY FROM DEBYE ENTROPIC GRAVITY 



Modified Newtonian dynamics (MOND) was proposed to explain the flat rotational curves of spiral galaxies. A 
great variety of observations indicate that the rotational velocity curves of all spiral galaxies tend to some constant 
value [l5j . Among them are the Oort discrepancy in the disk of Milky Way the velocity dispersions of dwarf 
Spheroidal galaxies [l?} and the flat rotation curves of spiral galaxies [l8j . These observations are in contradiction 
with the prediction of Newtonian theory because Newtonian theory predicts that objects that are far from the galaxy 
center have lower velocities. 

The most widely adopted way to resolve these difficulties is the dark matter hypothesis. It is assumed that all 
visible stars are surrounded by massive nonluminous matters. Another approach is the MOND theory which was 
suggested by M. Milgrom in 1983 [19]. This theory appears to be highly successful for explaining the observed 
anomalous rotational-velocity. In fact, the MOND theory is (empirical) modification of Newtonian dynamics through 
modification in the kinematical acceleration term 'a' (which is normally taken as a = v 2 jr ) as effective kinematic 
acceleration a c g = a/x(^), 

/ a n GM 

ao R 

where fj, = 1 for usual- values of accelerations and \x = ^-("C 1) if the acceleration 'a' is extremely low, lower than a 
critical value ao = 1Q -10 m/s 2 . At large distance, at the galaxy out skirt, the kinematical acceleration 'a' is extremely 
small, smaller than 10 -10 m/s 2 , i.e., a <C ao, hence the function /u(^) = Consequently, the velocity of star on 
circular orbit from the galaxy-center is constant and does not depend on the distance; the rotational-curve is flat, as 
it observed. 

Although MOND theory can explain the flat rotational curve, however its theoretical origin remains un-known. 
Thus, it is well motivated to establish a gravitational theory which can results MOND theory naturally. In this 
section, we are able to show that the MOND theory can be extracted completely from the Debye entropic gravity. 
This derivation further support the viability of Debye entropic gravity formalism. 

Again, we consider a spherical holographic screen with radius R as the boundary of the system. Combining Eqs. 
(3) and (4), and using the equivalence between mass and energy as well as relation A = AttR 2 , we obtain 

2tt , , GM , x 

—TD(x) = w . (29) 

Using the Unruh temperature formula ([8]). the above equation may be written as 

aD{x) = (30) 

Also, if we use the Unruh temperature formula in the definition of x, i.e. Eq. ([S]), and define ao as 

- 12Td 

a = — — , (31) 
irn 

then we obtain 

?r 2 ao 

X = IkT' ( 32 ) 
nt 7T 2 a GM 

aD(-) = (33) 

This is the MOND theory resulting from Debye entropic gravity. If we compare this equation with well-known Eq. 
(|28[). we see that we can define (x function as 

ti±) = D&°). (34) 
ao ba 

In what follows we show that this function satisfies the conditions similar to those of \x function in Eq. (|2"8")l . Let us 
examine Eq. (1331) in two limits of temperatures. First, we consider the limit corresponding to the temperatures large 
relative to the Debye temperature. In this case x <§C 1 (a 3> ao) we have D(x) = 1. Thus Eq. (|33|) reduces to 

a = -ftT- (35) 



Using the above result in Eq. (|30]) gives 
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Therefore, for strong gravitational fields, Eq. (|33|) turns into the standard Newtonian dynamics. As we discussed, for 
a > ao we have also A t (^-) = 1- We conclude that in the limit of a 3> a both D(x) and /j,(x) have the same behavior 
and become equal to 1. 

The second limit corresponds to the temperatures extremely smaller than the Debye temperature, T <C To, that 
is to say in the weak gravitational fields. In this limit, we have x ^> 1 (a <C ao), and the Debye function can be 
expanded as 

D(x) = - r -^—dy » — . (36) 
x J Q e y — 1 ox 



If we use the approximation (|36l) in Eq. (|33p , we obtain 



a 



(IQ 



GM 



« I - ) = #r- (37) 



Therefore, the Newtonian dynamics is modified for weak gravitational fields, e.g. at large distance from the galaxy 
center, namely at the galaxy out skirt. Thus the origin of the MOND theory can be understood completely in the 
framework of Debye entropic gravity. In this way we fill in the gap existing in the literature understanding the 
theoretical origin of MOND theory. 

IV. POISSON EQUATION FROM DEBYE ENTROPIC FORCE 

Finally, we obtain the modified Poisson equation by taking into account the Debye correction to the equipartition 
law of energy. We choose a holographic screen S corresponding to an equipotential surface with fixed Newtonian 
potential <po- We assume that the entire mass distribution given by p(x) is contained inside the volume enclosed 
by the screen and there are some test particles outside this volume. To identify the temperature of the holographic 
screen, we take a test particle and move it close to the screen and measure its local acceleration. The local acceleration 
is related to the Newton potential as 

a = -f(f>. (38) 
Substituting this relation into Unruh temperature formula, we get 

T _ ft|V0| 
2vr 

Using the above equation in the definition of x, we have 

T D 2wT D , , 

T h\V<j>\ 

Inserting (|39|) in Eq. (J7J , after using Eq. ([3]) for the number of bits on the holographic screen, we obtain 

M = i// W ^ dl - (41) 
Using the divergence theorem we can rewrite Eq. (|4ip as 

M = ^qJ V.[D{x)V<l>}dV. (42) 
On the other hand, for the mass distribution M inside the closed surface S, we have the relation 



(39) 



M= / p(x)dV. (43) 
Jv 

Equating Eqs. (|42]) and ([43]), we get 

f.[D(x)$(j>] = AixGp{x). (44) 
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This is the modified Poisson equation which is valid in all range of temperatures. For high temperatures, i.e. strong 
gravitational field (x <C 1) and hence D(x) = 1. In this case Eq. (PPfj) reduces to the standard Poisson equation, 

V 2 = 4irGp{x). (45) 

Thus, considering the gravitational system as a thermodynamical system and taking into account the Debye model 
for the modified equipartition law of energy, we see that not only Einstein equation and MOND theory but also the 
Poisson equation is modified accordingly. Clearly the modification of Poisson equation leads to modified Newton's 
law of gravitation. 



V. CONCLUSIONS 



In his work, Verlinde applied the equipartition law of energy as E = \NT on the holographic screen induced by the 
mass distribution of the system, and obtained the Einstein equations, Newton's law of gravitation and the Poisson 
equation. But we know from statistical mechanics that the equipartition law of energy does not hold at very low 
temperatures and it should be corrected. In this paper, we considered the Debye correction to the equipartition law of 
energy as E = ^NTD(x), where D{x) is the Debye function. Following Verlinde's strategy on the entropic origin of 
gravity, we obtained the modified form of the Einstein equations, MOND theory and the modified Poisson equation. 
Interestingly enough, we found that the origin of MOND theory can be understood from the Debye entropic gravity 
scenario. Since the MOND theory is an acceptable theory for explanation of the galaxy flat rotation curves, thus the 
studies on its theoretical origin is of great importance. This result is impressive and show that the approach here is 
powerful enough for deriving the modified gravitational field equations from Debye model. We also showed that in the 
temperatures extremely larger than the Debye temperature (very strong gravitational fields), the obtained modified 
equations turn into their respective well-known standard equations. The results obtained here further support the 
viability of Verlinde's formalism. 
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